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Abstract. 

By comparing radial velocities of radio bright compact H n regions with their H I absorption profiles, we discovered 
expanding shells of neutral hydrogen around them. These shells are revealed by absorption of the radio continuum 
emission from the HII regions at velocities indicating greater distances than the observed radial velocity. We believe 
that these shells are shock zones at the outer edge of the expanding ionized region. Additionally we found evidence 
for a velocity inversion inside the Perseus arm caused by a spiral shock, which results in a deep absorption line in 
the spectra of compact Hn regions behind it. 
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1. Introduction 

Models of the photodissociation of molecular clouds sur- 



rounding newly formed O and early B typ e stars (Hill & 
Hollcnbach, 1978| ; [Roger fc Dewdney, 199% show that the 
resulting H I will exist in a layered structure surround- 
ing the star. A dissociation front (DF) will first rapidly 



move through tire molecular cloud funning a broad Hi re- 
gion. As Lhe H n region expands a layer uf shucked H i will 

form just outside the ionization front (IF). Eventually the 
faster moving IF will catch up with the DF and all of the 
H I will be found in the shocked layer. The time when the 
IF and DF merge ranges from ~ 10 5 — 10 6 years depending 
upon the spectral type of the star and the densi ty of the 
surrounding molecular material (see Figure 7 of Roger & 
Dewdney, 1992). The broad Hi region has been detected 
in 21-cm line emission in a number of studies (e.g 



ber of cases we detect absorption features that are asso- 
ciated with shocked expanding layers of photodissociated 
gas. We focus on compact H n regions within the observed 
area of the Canadian Galactic Plane Survey (CGPS)Q 
(74.2° < I < 147.3°; -3.6° < b < 5.6°). The CGPS is 
part of the current international effort to map the Galactic 
plane in the 21-cm Hi line at ~ 1' resolution (Taylor 



1999; McClure-Griffiths ct al., 2001). At this resolution 



Roger. 



1982), but observational evidence for the shocked layer is 



currently limited to a study of the H II region Orion A by 



van der Werf fc Goss (1989) 



In this- 



paper we present a novel technique for the de - 
tection 1 and 3tudy of expanding II I structures around com- 

pact Hn regions that uses 21-cm Hi absorption spectra. 
The technique allows H I surrounding H n regions to be 
studied even if the H n region is not resolved. In a mini- 



mally Galactic H II regions will remain unresolved and the 
surrounding neutral material must be studied via absorp- 
tion. The technique presented in this paper, applied to 
such regions, will greatly increase the number of Hil re- 
gions studied using 21-cm absorption spectra. 

Hi absorption spectra towards radio bright Galactic 
Hn regions have been used primarily to probe the prop- 
erties o f the neutral ISM along the line of sight to the H n 
region ( Wcndker fc Wrigge, 1996 ; Normandeau, 1999 ). In 
these studies absorption features at velocities associated 
with the H n region itself are usually ignored. Another ap- 
plication has been in Galactic structure studies; Kuchar & 



Bania (1994) and Kuchar & Bania (1990) have shown how 



21-cm absorption spectra can resolve the distance ambi- 
guity towards inner galaxy H II regions. 



Send offprint requests to: R. Kothes 



1 CGPS data are publicly available via the Canadian 
Astronomy Data Centre (CADC; http://cadcwww. hia.nrc.ca) 



2 



Kothes & Kerton: Shocked H I shells around compact H II regions 



The H I absorption spectrum towards a galactic H II re- 
gion can also be used to study the kinematics and struc- 
ture of neutral material associated with the Hll region. 
Such absorption studies have tended to be limited to a few 



very bright Hll regions such as DR 21 (Roberts, Dickel, & 
Goss, j997l), Orion (|van der Werf fc Goss, 1989|) and W3 



van der Werf & Goss, 1990). More recently Lebron et al 
(200l|)Jused Hi absorption as part of a multiwavelength 
study of the blister Hll region Glll.61+0.37. 

In §2 we describe the data used and outline the data 
processing techniques. Particularly we focus on a new spa- 
tial structure filtering technique. In §3 we show the ab- 
sorption spectra of all our Hll regions. A general model 
for the structure of an absorption spectrum is presented 
and applied to our observations in §4 followed by a general 
discussion in §5. Finally, a summary is presented in §6. 



2. Observations and Data Analysis 

2.1. H I Line Data from the CGPS 

Hi 21-cm line data were obtained using the DRAO syn- 
thesis telescope (Landecker et al., 2000) as part of the 
Canadian Galactic Plane Survey (CGPS) ( [Taylor, 199i| . 
The Hi data reduction techniques used in the CGPS are 
described in detail in Higgs (20"00| ). The CGPS consists of 
data cubes with a size of 5.12° x 5.12°. For further analysis 
we extracted 1.4° x 1.4° fields for each Hll region studied. 



Name 


T c e 


vlsb 


Reference 


other names 




[K] 


[km/s] 






G75.77+0.35 


450±23 


-8.5 


(1) 




G75.83+0.40 


750±38 


-4.1 


(1) 


Sh 105 


G76.15-0.29 


100±8 


-28.2 


(1) 




G76.19+0.10 


55±6 


-2.1 


(2) 


IRAS 20220+3728 


G80.35+0.73 


100+11 


-62.4 


(1) 




G80.94-0.13 


374±21 


-2.1 


(1) 


IRAS 20375+4109 


G83.94+0.77 


143±9 


-8.5 


(1) 




G85.25+0.02 


69±7 


-31.2 


(1) 




G85.69+2.03 


41±5 


-75.0 


(4) 


IRAS 20446+4613 


G90.24+1.72 


32±5 


-60.9 


(3) 


Sh 121 


G92.69+3.08 


48±6 


-6.4 


(2) 


IRAS 21078+5211 


G93.86+1.56 


44±5 


-59.5 


(2) 


IRAS 21202+5157 


G96.29+2.60 


52±5 


-94.7 


(3) 


Sh 127 


G97.51+3.17 


72±5 


-71.1 


(2) 


Sh 128 


G105.62+0.34 


71±5 


-52.1 


(2) 


Sh 138 


G108.20+0.58 


259±14 


-55.9 


(1) 


Sh 146 


G108.37-l.06 


58±5 


-53.1 


(3) 


Sh 149 


G108.76-0.95 


150±9 


-49.1 


(1) 


Sh 152 


G110.11+0.05 


225±12 


-51.7 


(1) 


Sh 156 


Glll.62+0.37 


87±6 


-63.4 


(1) 


Sh 159 


G112.23+0.24 


82±7 


-41.8 


(1) 


Sh 162 


G114.02-l.45 


63±5 


-49.0 


(5) 


IRAS 23365+5953 


G115.79-l.58 


33±4 


-40.8 


(1) 


Sh 168 


G132.16-0.73 


79±5 


-56.6 


(1) 


IRAS 02044+6031 


G137.75+3.80 


77±5 


-47.0 


(5) 


IRAS 03030+6229 


G138.49+1.64 


61±5 


-34.6 


(1) 


Sh 201 



2.2. Data Analysis 



Table 1. Characteristics of the compact Hn regions in 
our sample. Tjf represents the peak continuum bright- 
ness temperature in the CGPS 1420 MHz data and vlsr 
the radial velocity. References for the radial velocities: (1) 
Lockman (19891), ( 2 ) |Bronfman et al. (1996[), (3) |Brand & 



To dctl crminc H I emission and absorption profiles we Blitz (1993|), 7 1) |Woutcrloot fc Brand (1989| ), (5) 



have to separate the Galactic background emission from 
the small scale structures around and on the compact 
Hn regions we plan to study. For that purpose w e used 
the "background filtering technique" invented by Sofue 
fc Reich (1979| ). This method was developed for separat 
ing smooth Galactic background emission from small scale 
emission structures. We have modified the technique to be 
useful for negative features as well as for positive ones. 

For each velocity channel we used the following proce- 
dure: 



Brunt & 

Kerton (2002| ). Ref. (1) contains recombination line veloc 



Tg convolution to a beam 8 



Tb — Tg — > ST B 
T b x ST B -> UT B 

T b : UT B < 
5T B : UT B > 

Tlarge Convolution to 

Tb T[ ar g e ► T sma n 



T 1 - 

1 B — 
1 B 



Here Tb is the original brightness temperature distri- 
bution with subtracted continuum emission, Ti arge rep- 
resents the large scale background emission, and T sma u 
the small scale structures including the absorption of the 
compact Hn regions. 

For our Hi data this technique creates one Hi data 
cube representing large scale structures which provides in- 
formation about the emission of the smooth background 



itics and all other references radial velocities of molecular 
lines. 



component for each source. We also get one H I data cube 
containing small scale structures. Absorption profiles to- 
wards the compact H n regions and emission profiles of 
possible H I structures around them are derived from this 
cube. 

3. Absorption profiles of compact H n regions 

3.1. The source sample 

We obtained absorption profiles towards all Hn regions 
in the CGPS survey area which are unresolved in the 
1420 MHz continuum data and have a flux density greater 
than 200 mJy (T(, ~ 30 K) therein. We obtained radial ve- 
locities for these H n regions from various sources as listed 
in Table [l]. Sources with positive radial velocities were re- 
moved from the sample because of the distance ambiguity 
in the inner Galaxy. All sources in our sample are bright 
enough in the CGPS to show significant absorption by 
foreground neutral hydrogen. 
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3.2. Obtained absorption profiles 

We used the background filtering technique described in 
Sect. 2.2 to separate the broad structures from the com- 
pact absorption features. We used a beam O of 5' for all 
sources and checked all final data cubes carefully by eye. 
For three stronger sources, G75. 77, G75.83, and G80.94, 
there were still some faint features from the absorption 
visible in the background component. For these we used a 
beam of 10' which then successfully separated small scale 
and large scale structures. We created the absorption spec- 
tra in Figs. [| to U by simply plotting the value at the peak 
position of the radio continuum source in each H I chan- 
nel of the H I cubes as a function of the radial velocity. 
One naively expects to see absorption of the H II region's 



continuum emission by the nentral hydrogen in Lhe fore- 
ground up Lo Us radial velocity maybe wiLh a Gaussian 

slope with a width of a few km/s towards higher nega- 
tive velocities, assuming a flat rotation model. In contrast 
we see at least one additional absorption line well beyond 
the radial velocity in all but three of the 26 observed H II 
regions. The additional absorption lines are at higher neg- 
ative velocity than the radial velocity of the source itself 



The shocked layer of H I initially expands with the H n 
region at velocities ^10 km/s. For extremely young sys- 
tems, where there has been little or no expansion of the 
Hn region, the shocked layer of Hi will not be present 
(e.g. Dewdney & Roger, 1982). By the time IF merges 
with the more slowly moving DF the expansion velocity 
of the shocked layer will have dropped to ~ 3 — 4 km/s. 
The merging of the two fronts occurs on a time scale of 
10 5 — 10 6 years depending upon the density of the molecu- 
lar material and the spectral type of the exciting star. The 
unshockcd H I layer is not expected to have any systematic 
velocity shift relative to the molecular material if it is fully 
confined by the molecular material. Observations that do 
show expansion of the unshocked dissociated H i are all 
associated with blis ter Hil regions ( Roger, 1982] ; van dci 



Wcrf fc Goss, 1989). The observed expansion velocities of 



2-6 km/s in these cases are most likely due to the expan- 
sion of the heated dissociated gas into the surrounding low 
density medium. 

Gas temperatures are expected to be < 100 K within 
the unshocked material; newly shocked material may be 
hotter but is expected to cool rapidly to T= few x 100 K 



indicat ing material which is moving 



towards us relative because of effici ent cooling in the shocked layer ( [Roger fc 



to the I H II region , Considering the fact that we find this 

feature in almost all of the observed compact H n regions, 
we can assume that it is a symmetric expanding shell of 
neutral hydrogen. Since we are dealing with compact H II 
regions which are unresolved in our data we do not find 
those expanding shells around them in emission because 
the e mission is buried together with the foreground ab- 



Dcwdney, 1992| ). Observational support for these values 



2 mis 
tioL 



come from ISO observations of H2 lines towards photodis 
sociation regions around S 140 (Timmcrmann ct al., 1996) 
and S 106 (van den Anckcr, Tielens, & Wesselius, 2000) 
which indicate T~ 500 K, and from Hi 21-cm line emis- 
sion observations of a compact H 11 region in the nebular 
complex GGD 12-15 which indicate T~ 300 K flGomcz et 



sorption within the beam. 



al, 1998) 



We should note at this point that the additional ab- 
sorption line could also be caused by a double valued ro- 
tation curve. Our sample was chosen to avoid the velocity 
ambiguity in the inner Galaxy, but there is also evidence 
for a velocity inversion in the Perseus arm caused by a 



spiral shock (Roberts, 1972). However, since we observe 
this phenomenon for almost all of the compact Hil re- 
gions independent of the location and distance we assume 
this to be a minor factor. Nevertheless, we will discuss this 
in more detail in Sect. 5. 



4. Theory and Modelling 

4.1. Dissociated molecular gas around massive young 
stars 

In this section we summarize some of the basic results 



The column density of shocked H 1 is expected to grow 
as the system evolves, and has t ypical values of 10 20 — 10 2 1 
over the lifetime of the system ( Roger fe Dewdney, 1992 ). 
The unshocked layer is expected to reach a peak column 
density of ~ 1 — 3 x 10 21 with little dependence on the 
spectral type of the exciting star and the density of the 
molecular material ( Roger fc Dewdney, 1992 ). 



4.2. A simple model for expanding H 1 shells 

To calculate observational properties of the expanding H 1 
shells around the compact H 11 regions we have developed 
a simple model based upo n Hill fc Hollcnbach (1975 ) and 
Roger fc Dewdney (1992] ). An early type star, located in 



of the models of Hill & Hollcnbach (1978) and Roger & 
Dewdney (1992) regarding the morphology, velocity, tem- 
perature and column density structure of the dissociated 
gas surrounding a newly formed OB star. Essentially these 
models demonstrate that the dissociated gas will have a 
layered structure consisting of a broad unshocked layer 
formed by a dissociation front (DF) and a thinner shocked 
layer of Hi located just outside of the expanding ioniza- 
tion front (IF). 



the centre of a molecular cloud, creates an expanding H 11 
region inside an area of dissociated molecular material (sec 
Fig. ^|) . At the interface we have a narrow zone of shocked 
neutral material expanding together with the ionization 
front. For simplicity we assume the whole structure to 
be symmetric. The object is located at a distance d from 
the Sun which translates to a radial velocity Vd assum- 
ing a flat rotation curve with v@ = 220 km/s and Rq — 
8.5 kpc. We further assume that the object is located in 
the outer Galaxy so that the velocity Vd is negative and 
a higher negative velocity implies a greater distance from 
the Sun. The part of the shocked H 1 shell moving towards 
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us would be shifted in its radial velocity to Vd — v exp , with 
v exp representing the expansion velocity of the Hi shell, 
while the velocity of the receding shell would be shifted to 
v d + v ex P (see Fig. ^|). In this case the receding Hi, even 
though it is located behind the Hn region, would be ob- 
served among Galactic H I located between the object and 
the Sun. The approaching part of the shell, even though 
located in front of the H II region, would appear among 
the material which is located behind the object in our H I 
data (see Fig. ^) . 

Assuming a constant spin temperature Ts within the 
shocked neutral material the observed brightness temper- 
ature Tb at the position of the H n region would be: 



(1) 



(2) 



T B (v) = TJ(1 - e- T W) 

for the receding part and: 

T B {v) = TJ(1 - e- T{v) ) + T c (e~ T(v) - 1) 

for the approaching part, where the second term repre- 
sents the absorption of the H n region's continuum emis- 
sion by the approaching Hi shell. Here T c is the Hn re- 
gion's effective continuum brightness temperature in our 
1420 MHz data and r is the optical depth within the 
H I shell. Tg is the effective spin temperature. Since these 
compact H II regions are unresolved in our data the peak 
brightness temperature is not the real brightness temper- 
ature of the source. Tq is Tq multiplied by the ratio of the 
actual source area to the observed source area, that is: 



Pi 2 



e 



where 8^ is the HPBW of the observation and 0s the 
angular diameter of the source. For a source with a di- 
ameter of 10" and a brightness temperature of 1000 K for 
example we would measure a peak of only 30 K. The same 
effect applies to the spin temperature. 

If we now create H I absorption line spectra as de- 
scribed in Sect. 3.2 we will find continuum absorption of 
the H II regions radio emission by the foreground material 
up to its radial velocity. Since we find the receding part 
of the expanding H I shell at a radial velocity of Vd ± v exp 
it will appear as an emission line superimposed on the 
absorption structure of the foreground material making 
it hard to detect. The approaching part, however, has a 
radial velocity of Vd — v exp . We expect to find its related 
line separated from the foreground absorption at a higher 
negative velocity. 

Since our sample consists of sources which are signif- 
icantly smaller than the beam of the observations, we 
can only speculate about the structure of neutral mate- 
rial around these compact H n regions. We can resolve 
different layers with different expansion velocities but we 
do not get any information about the spatial distribution. 
Our method favours the part of the expanding shells which 
sees most of the H n region's radio continuum emission. In 
this case small high density knots are smoothed out and 
thus not detectable. 



Name 


[km/s 




T 


T 

T s /T c = 0.2 


HPBW 

[km/s] 


G75.77+0.35 


19.2±0 


6 


0.17 ±0.01 


0.21 


1.7 


G75.83+0.40 


25.6±0 


7 


0.05 ±0.01 


0.07 


1.7 


G76.15-0.29 


7.7±3 


2 


0.45±g-gg 


0.61 


4.9 


G76.19+0.10 


5.1±0 


8 




1.15 


2.8 


G80.35+0.73 


5.4±2 


9 


0-64+jjif 


0.89 


5.1 


G80.94-0.13 


2.7? 






_ 


_ 


G83.94+0.77 


11.1±2 


2 


0-28tn 3 2 


0.37 


5.6 


G85.25+0.02 


10.8±2 


4 


0-49±H 8 6 


0.66 


2.4 


G85.25+0.02 


5.9±2 


4 


0.23 ± 0.03 


0.30 


1.2 


G85.69+2.03 


2.5±0 


5 


0.26 + °iS 


0.34 


3.8 


G90.24+1.72 


23.0±0 


5 


0.9lt°j 3 


1.37 


2.4 


G92.69+3.08 


0.8±0 


8 


> 2.00 


> 2.00 


4.2 


G93.86+1.56 


11.0±0 


8 


0.821°^ 


1.21 


6.8 


G93.86+1.56 


2.2±0 


8 




0.62 


2.9 


G96.29+2.60 


_ 






_ 


_ 


G97.51+3.17 


9.4±0 


8 


0.70+°'£° 


1.00 


4.2 


G97.51+3.17 


2.2±0 


8 


1.60l°;« 


>2.00 


3.5 


G105.62+0.34 


4.0±0 


8 


0.28 ±0.03 


0.36 


1.2 


G105.62+0.34 


1.5±0 


8 


0-681°;™ 


0.97 


4.2 


G108.20+0.58 


_ 






_ 


_ 


G108.37-l.06 


0.6±1 


3 


> 2.00 


> 2.00 


3.6 


G108.76-0.95 


4.3±0 


9 


l-76±g-« 


> 2.00 


5.8 


G110.11+0.05 


6.2±1 


2 


0-78±° ^ 


1.13 


3.2 


Glll.62+0.37 


_ 






_ 


_ 


G112.23+0.24 


9.4±1 


5 


0.16 ±0.02 


0.21 


1.0 


G114.02-l.45 


5.4±0 


5 


o.89ig;2 


1.33 


1.7 


G115.79-l.58 


4.0±1 


5 


1.46ig;» 


> 2.00 


1.6 


G132.16-0.73 


2.5±4 


8 


1 Cj4+0-35 
1M -0.26 


> 2.00 


3.7 


G137.75+3.80 


6.1±0 


5 


U - Oi -0.05 


0.70 


1.3 


G138.49+1.64 


4.3±1 





> 2.00 


> 2.00 


3.1 



Table 2. Derived parameters for the compact Hn regions 
in our sample. Sources with two additional absorption 
lines are noted twice. 



4.3. Application to the data 

We can now calculate characteristics of the expanding H I 
shells by fitting Gaussians to the individual absorption 
spectra. Four examples are shown in Fig. ^. A compari- 
son of the H II region's radial velocity with the position of 
the additional absorption line gives the expansion veloc- 
ity of the H I shell. By transforming eqn. || we obtain an 
expression for the opacity r within the expanding shell: 



-In 



T C -TI 



(3) 



We do not have any information about the spin tem- 
perature within the expanding H I shell since its emission 
is superimposed on the absorption profile. Assuming that 
the spin temperature is significantly lower than the con- 
tinuum brightness temperature of the H n region we can 
neglect it in eqn. || and calculate a value for r. We should 
note at this point that the effect of Ts on the resulting 
value of r increases with high Tb to Tq ratios. As an ex- 
ample we also calculated the resulting r when the spin 
temperature is 20 % of the continuum temperature. The 
results for the individual Hil regions are listed in Table 



Kothes & Kerton: Shocked H I shells around compact H II regions 



5 



g. For r > 2.0 we have indicated a lower limit since the 

quantity ^T^™ is smaller than our uncertainties in this 

1 c 1 s 

case. Three of the sources do not have additional absorp- 
tion lines at higher negative velocities and we could not 
obtain a satisfactory multiple Gaussian fit to the absorp- 
tion profile of G80.94 due to its complex structure. 



5. Discussion 

5.1. The nature of the additional absorption line(s) 

The discussion about the evolution of the layer of shocked 
neutral hydrogen expanding with the ionizing shock front 
around young stellar objects in Sect. 4.1 gives us con- 
straints on the characteristics of these features. We find a 
maximum expansion velocity of about 10 km/s. An upper 
limit of about 1.0 for the opacity r can be deduced from 
the maximum Hi column density, Nh = 10 21 cm -2 , and 
the typical spin temperature of a few 100 K by using the 
relation: 



iV, 



Hi 



with 



1.823 • 10 18 • To 



r(v)dv 



(4) 



~(v)dv = 1.06 ■ t ■ 5v 



for a Gaussian distribution. Here Sv is the HPBW of the 
absorption line, which usually lies around a few km/s and 
r is the opacity at the peak of the line as calculated in 
Table |^. For an examination of those constraints within 
our data we plotted the opacity r as a function of the 
expansion velocity for all available sources in Fig. ||. For 
those sources where we found double structures we plotted 
both results to allow for multiple lines of different origin. 

pairs 



u exp 



We found three well separated groups of r 
in the diagram. Most sources are located in the lower left 
corner. Their t and v exp values agree very well with the 
constraints we laid out for the expanding layers of shocked 
Hi. There seems to be a sharp upper limit for the expan- 
sion velocity at about 11 km/s and an upper limit for r 
of about 1.0. 

There is another group of sources with low v exp and 
very high r. We should note at this point again that the 
t values in this diagram were calculated by assuming that 
Ts << Tc- The effect of T$ on the final r is bigger for r 
values which are already high by assuming T$ « Tc (see 
Table^) . Hence the sources located in the lower left corner 
are not that much affected by higher T$ values than those 
which are already high above this area. This means that 
we cannot explain the sources located high above r = 1 
with expanding layers of shocked Hi. 

We believe that the absorption features we are detect- 
ing in the second group indicate a possible distance am- 



biguity in the Perseus spiral arm. According to Roberts 



(1972) the spiral shock causes the radial velocity to drop 
by 20 - 30 km/s at the beginning of the Perseus arm and 
slowly rejoins the flat rotation curve after about 1 kpc. 



This creates a velocity inversion inside the Perseus arm. 
Thus inside the arm the distance is increasing with in- 
creasing radial velocity as opposed to the flat rotation 
curve. According to Roberts (1972) the density in the 
closer part of the Perseus arm is significantly enhanced by 
the spiral shock which would explain the high r values we 
obtained. For this interpretation we would expect a strong 
dependence of the radial velocity of this absorption line on 
Galactic longitude. Because of the changing viewing angle 
to the Perseus arm the radial velocity should increase with 
increasing Galactic longitude. To examine this dependence 
we plotted the vlsr at the center of the absorption line 
as a function of Galactic longitude, for the sources with 
high opacities in Fig. || (filled stars). We added the last ab- 
sorption lines from the spectra of the three sources which 
have no additional line at higher negative velocities and 
the source G90.23, one of the high velocity sources in the 
right part of Fig. || (open stars). We believe this source 
belongs to this group as well, because for a high velocity 
wind the r is too big, and for the dissociated material, 
shocked or unshocked, the expansion velocity would be 
too high. 

Apparently there is a strong relation between the ra- 
dial velocities and the Galactic longitude. There are only 
two sources which do not follow this relation. One of those 
exceptions is the source G92.68 located close to the upper 
left corner in Fig. || which is of course local based on the 
low negative radial velocity and should not show the spiral 
shock in the Perseus arm at all. It has one of the lowest 
expansion velocities (0.8 km/s) which is comparable to the 
uncertainty in our data. So this strong absorption line is 
most likely due to a dense H 1 cloud in the foreground. The 
second exception is the source G132.16 in the lower right 
area of Fig. |9[ The r is too high to be explained by the ex- 
panding shocked neutral hydrogen layer, but as mentioned 
in Sect. 4.1, there are rare occasions where the shell of dis- 
sociated material outside the shock front has a systematic 
velocity of 2-6 km/s for blister Hn regions. These shells 
have higher H 1 column densities and lower spin temper- 
atures so we would expect higher opacities. Additionally 
G132.16 shows another strong absorption line in its spec- 
trum at about -40 km/s which would fit nicely to the spiral 
shock model. 

There are 4 sources in our sample which show two ad- 
ditional absorption lines in their spectra instead of one. 
G97.51 shows one line with the characteristics of the spi- 
ral shock and the expansion velocity of the shocked neutral 
hydrogen is apparently high enough so that both absorp- 
tion lines are well separated. The sources G85.25, G93.86, 
and G105.62 all have one high and one low velocity com- 
ponent both of which have rather low opacities. We believe 
that the higher velocity component is the result of absorp- 
tion by the shocked neutral hydrogen while the component 
with the lower velocity was produced by the unshocked hy- 
drogen. This would indicate that all three are blister H 11 
regions. 

Finally, two of the observed H11 regions, G75.77+0.35 
and G75.83+0.40, have a very high v exp and low r. Since 
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both of these sources are very close together on the sky 
(only 5' apart), one possibility is that we are seeing unasso- 
ciated blueshifted foreground Galactic gas. We think this 
is unlikely though since the observed absorption line is 
seen at different velocities in each case (v ra d = —27.7 km/s 
for G75.77 and v rad = -29.7 km/s for G75.83). Instead 
we believe we are observing a fast neutral outflow from 
these objects, perhaps formed from a recombining ion- 
ized wind. These two objects are by far the most ener- 
getic Hn regions in our sample, and previous 12 CO ob- 
servations have detected high velocity molecular gas, with 
v ~ 10 — 50 km/s, towards these objects ( [Shepherd fc 



Churcl well, 1996). Models of magnetic fields in Hn re- 
gions suggest that for certain magnetic field alignments 
fast neutral winds (v > 20 km/s) can form from recom- 
bining ionized flows that have broken out of the molecular 



cloud ( Franco, Tcnorio-Taglc, fc Bodenheimcr, 1989| ). In 



addition, a fast neutral outflow with v < 100 km/s has 
been observed in the high-mass star forming region DR 21 
( Russell et ah, 1992 ). From this, we conclude that we are 
seeing in absorption the blueshifted part of a fast neutral 
outflow coming from these objects although the linewidth 
is surprisingly narrow in both cases. The velocity shift in 
this case is the outflow velocity projected along the line of 
sight. 

5.2. Where is the receding shell? 

We identified most of the additional absorption lines as 
the absorption of the H n region's radio continuum emis- 
sion by the expanding layer of shocked neutral hydrogen 
around the ionized gas. Now we have to discuss the re- 
ceding parts of these symmetric layers. By utilizing eqns. 
(1) and (2) we can calculate the amplitude ratio of the 
absorption line of the approaching shell to the emission 
line of the receding shell, which results in 



T B (v d - v ea 



T B {v d + v exp ) 



\T S -T C \ 
To 



This ratio would be 1 if the radio continuum temperature 
is twice the spin temperature. For the compact H n regions 
we should expect radio continuum brightness tempera- 
tures of up to 1000 K, maybe even 2000 K, at 1420 MHz, 
while the spin temperature should be around a few 100 K 
up to a maximum of 1000 K. Since most sources show 
rather complex foreground absorption spectra it is very 
difficult to decide if a structure is an emission line su- 
perposed on the foreground absorption or just the lack of 
absorption at this position. In four cases, G76.15, G76.19, 
G93.86, and G114.07, we may see the receding shell in 
emission. All of these sources show a possible line at about 



the expected radial velocity va 



We get amplitude 



ratios between the absorption line of the approaching shell 
and the emission line of the receding shell of about 3:1 for 
G76.15 and G114.07 and about 3:2 for G76.19 and G93.86. 
This would translate to T c to T s ratios of 4 and 2.5 respec- 
tively. 



For those objects where we do not see the emission line 
one can calculate the T c to T s ratio required to make the 
line undetectable. The highest T c to T s ratios we derive 
are for G80.35 and G85.25. These sources show neither 
emission nor absorption at the expected radial velocity 
for the receding shell and we need ratios of around 6 so 
that the amplitude of the receding shell is hidden in the 
noise. This is not a problem for the typical values for T c 
and T s we discussed above. 

5.3. G76. 15-0.29, a special case? 

For none of the sources in our sample we would expect 
to observe either of the neutral layers around them in 
emission because these sources are unresolved and thus 
this emission should be buried within the absorption lines. 
Especially the thin shocked layer should be impossible to 
detect in our data. However, around the source G76.15, 
which is slightly extended in our data, we found an incom- 
plete shell of Hi emission centered on its recombination 
line velocity (see Fig.p^|). This Hi shell has a u-shaped 
structure open to the north-east. In the north-west and 
the south-east it seems to be located right outside the H n 
region while there is a gap between them in the south-west. 
We derived emission profiles for this structures at differ- 
ent radii from the center of the Hn region going to the 
north-west (see Fig.[ll]) . The separation of the two emis- 
sion peaks indicate an expansion velocity of 5 - 6 km/s. 
The existence of the emission gap between the H n region 
and the H I shell in the south-west demonstrates that the 
Hi shell cannot be the layer of shocked neutral hydrogen, 
since this layer should be located right outside the ioniza- 
tion front. Hence, this expanding shell must be the layer of 
unshocked dissociated molecular material. The only expla- 
nation for the systematic velocity difference of 5 - 6 km/s 
is that the Hn region is not completely surrounded by 
molecular material. Thus we can expect flows of neutral 
hydrogen towards low density areas giving the material 
a velocity along the line of sight. G76.15 must also be a 
blister H II region. 

6. Summary 

We developed a background filtering technique that allows 
the detection of small-scale positive and negative features 
within a data set. The filtering was applied to CGPS Hi 
observations towards 26 compact Hn regions in order to 
obtain absorption and emission spectra. 

In all but three of the spectra we observe absorption 
lines at more negative velocities than would be expected 
from the absorption of continuum emission by foreground 
material. A simple model of an expanding H n region with 
surrounding shocked and unshocked H I zones was devel- 
oped to interpret the observations. For the majority of the 
observations, the physical properties (expansion velocity, 
column density) of the shells derived from the absorption 
spectra are consistent with them being caused by shocked 
expanding shells of H I. 
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Some of the additional lines we observe are caused by 
a large scale Galactic shock associated with the Perseus 
arm of the Galaxy. Evidence for this comes from the way 
in which the velocity of the observed absorption features 
vary with longitude. 

Finally, two of the observed H n regions show evidence 
for high velocity neutral flows such as the flows observed 
in emission towards DR 21. 

This study has shown how absorption spectra can be 
used to examine neutral material surrounding H n regions 
even if the Hn regions are unresolved. The techniques 
that have been developed should prove to be useful as 
radio observations of the Galactic plane at 1' resolution are 
extended beyond the CGPS region to cover the majority 
of the plane. 
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Fig. 1. Hi absorption spectra of six compact Hn regions. These spectra were taken at the peak position of the 
continuum source. Radial velocities are indicated by a dashed line. Upper lines show the emission profile at the 
position of the Hn regions taken from Ti arge and the lower line the absorption profile taken from T sma u. Arrows 
indicate the absorption features of interest. 

Fig. 2. see Fig. 

Fig. 3. see Fig. [j], positive structures in the absorption profiles of G96.29 and G97.51 are filaments coincidentally at 
this position and not related to the H n regions 

Fig. 4. see Fig. [j], positive structures in the absorption profiles of G112.23 and G115.79 are filaments coincidentally 
at this position and not related to the H n regions 

Fig. 5. see Fig. |l| 

Fig. 6. Left: A sketch of a compact Hn region (H + ) with an expanding layer of shocked neutral hydrogen (H*) as 
described in Sect. 4.1. The dissociation front (DF) and ionization front (IF) are indicated, as is the layer of the 
unshocked Hi (H). Right: A diagram of the radial velocity as a function of the distance from the Sun at a longitude 
of 76° for a flat rotation model with vq — 220 km/s and Rq — 8.5 kpc. Velocity shifts of an expanding shell with 
v e xp = 10 km/s for an H n region with a radial velocity of Vd = —40 km/s are indicated. 

Fig. 7. Four example absorption profiles with the fitted Gaussians. The histograms represent the observed brightness 
temperature distribution and the solid line the fitted Gaussians. 

Fig. 8. Diagram with calulated minimum r values for all H n regions as a function of the expansion velocity. Three 
groups of sources with possible different origin for the additional absorption line are marked. 

Fig. 9. Diagram of the radial velocity as a function of Galactic longitude for sources with high calculated r values. 
Filled stars represent sources which are located in the box marked "spiral shocks" in Fig. ||. Open stars represent the 
last absorption line of those sources which have no additional absorption line at a velocity v < vlsr an d the source 
G90.24. 

Fig. 10. The Hi shell found in emission around the Hn region G76.15 integrated between -22 and -30 km/s. 
Fig. 11. Emission profiles of the Hi shell around G76.15 taken at different radii. 



This figure "fl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "f2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "f3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "f4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "f5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "f6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "f7.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "f8.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "f9.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "flO.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



This figure "fll.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0205032vl 



